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Abstract. We consider a simple market where a vendor offers multiple variants of a certain product and 
preferences of both the vendor and potential buyers are heterogeneous and possibly even antagonistic. 
Optimization of the joint benefit of the vendor and the buyers turns the toy market into a combinatorial 
matching problem. We compare the optimal solutions found with and without a matchmaker, examine the 
resulting inequality between the market participants, and study the impact of correlations on the system. 
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CT 1 1 Introduction 



The study of complex economic systems has attracted at- 
tention of many physicists. They have contributed to the 
field with several highly simplified, yet influential, mod- 
els such as the minority game percolation ||, scaling 
models of financial markets ip§[, as well as with a set of 
useful tools and insights [J, [5|, la Lfl] ■ 

Adopting a simplifying point of view characteristic for 
the works mentioned above, in this paper we focus on the 
interactions between consumers and producers. These in- 
teractions represent a classical example of decision-making 
under uncertainty [9( where the limited information avail- 
able to the contract participants results in a risk of mak- 
ing a wrong decision. In the standard economic litera- 
ture, problems related to the interactions of consumers 
and producers are as diverse as the research of consumer 
behavior the question of trust [10], the economics of 
information [TTI. [l3| ■ and the behavior of entire firms and 



industries lJTjjT 

This work is particularly motivated by the classical 
stable-marriage problem [IH, in which N men and N 
women all have their individual preferences and are to 
be matched one-to-one. Almost inevitably, it's impossi- 
ble to satisfy everyone and hence stable matchings (where 
no one has the possibility to exchange the assigned part- 
ner for a better one) or the optimal matching (where the 
total satisfaction is maximal) are of interest. The stable 
marriage problem has implications in many economic and 
social systems. It can represent matching of job seekers 
and employers or that of lodgers and landlords; it is also 
a metaphor for problems in logistics [l7| and in online 
marketing [181 ]. 

We study a situation where a certain product is avail- 
able in multiple variants and the preferences of both the 



buyer and the vendor for each of the variants can be rep- 
resented by numbers (the higher the number, the more 
appreciated the variant). The matching of the buyer and 
the vendor is then achieved by the selection of a single 
variant to be delivered. In the given framework, we first 
study outcomes achieved with the help of an external 
matchmaker — an idealized agent supervising the market 
and having perfect information about all the preferences. 
In particular, we investigate the inequality between profits 
enjoyed by the two involved parties and how correlations 
of the preferences influence the system's behavior. 

For comparison, we study two simple matchmaker-free 
models of variant selection. In the first matchmaker-free 
model, the vendor makes consecutive offers and the buyer 
decides whether to accept an offer or not. Our results show 
that while this approach results in a small decrease of 
the total satisfaction, it considerably decreases the buyer- 
vendor inequality. In the second matchmaker-free model, 
the buyer is searching for the optimal variant by himself. 
We study the optimal number of examined variants and 
show that under some conditions, this number may be 
infinite: the buyer is tempted to search forever. Numer- 
ical simulations of the studied models are in most cases 
accompanied with approximate analytical results. 



2 Trading under matchmaker's supervision 

We assume that a given product is available in N differ- 
ent variants which can be prepared by a vendor and fulfill, 
to a greater or lesser degree, needs of a given buyer. The 
buyer's utility from purchasing variant a is denoted by 
x a and the vendor's utility from providing this variant is 
denoted by y a . The matchmaker optimizes the joint ben- 
efit by maximizing the total utility u a {x a , y a ). Obviously, 
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the system's behavior depends on the choice of the utility 
function and on the nature of the utilities x a , y a — we shall 
study different settings in the following sections. 



2.1 Linear utility function 

The simplest form of the total utility is 

u a (x a , y a ) — -\- y a ( 1 ) 

where both utilities are merely summed with equal weights. 
In addition, we assume that both x a and y a are random 
variables drawn from the uniform distribution in the range 
[— 1, 1], U{— 1, 1), and that they are uncorrelated. The dis- 
tribution f(u a ) then has the tent-shaped form 



/(««) 



(2 
(2 



U a )/i 



[-2;0), 
[0;2]. 



The probability that a randomly selected variant has the 
total utility greater than u a is P(u a ) := f f(u')Au'. 
Since utilities of different variants are mutually indepen- 
dent, the largest utility u m :— max* =1 u a has the distri- 
bution 

g(u m ) = Nf(u m )[l - Piu^f- 1 . (3) 

Here the factor N appears because any of N variants can 
have the largest utility, the factor f(u m ) is the occurrence 
probability of u m , and the factor [1 — P(u m )] JV_1 is the 
probability that the remaining N — 1 variants have utili- 
ties lower than u m ; Eq. ([3]) is also known as the extreme 
statistics of the random variable u m . We study the model 
by computing (u m ) (by (x) we denote the average of x 
over all possible realizations). Since P(u a < 0) = 1/2 and 
u m < only when all N variants have u a < 0, it follows 
that P(u m < 0) = 2~ N . Hence, assuming N ^> 1, we 
can confine our computation to u m > where f(u m ) = 
(2 - u m )/A and P{u m ) = (2 - u m ) 2 /8. When N is large, 
P((u r n)) is small and thus we can use the approximation 
1 — P(u m ) w exp[— P(u m )]. After replacing the lower in- 
tegration bound in (u rn ) with — oo (this is again justified 
by the negligible probability of u m < 0) we obtain 



2N / 2ttN 2 
N~^l ~ V (N- l) 3 



V N 



(4) 



where we neglected terms of order 0(1/ N) and higher. 
We see that as N increases, (u m ) rapidly approaches its 
upper bound — the difference scales with TV -1 / 2 . Numeri- 
cal computation of (u m ) shows that the relative error of 
Eq. (|4]) decreases fast with N: it is less than 1% already 
for N = 17. 

Apart from the optimal total utility u m , the inequality 
between the vendor and the buyer is also of interest. If 
variant (3 maximizes u a , we say that the buyer- vendor 
inequality is \xp — yp\ and denote its expected value by A. 
To compute A, one needs to realize that for any given u a , 
the term \x a — y a \ ranges from to 2 — u a (this maximal 
difference is achieved when one of the two utilities is 1 



and the other is u a — 1). Since x a and y a are uniformly 
distributed, all possible values of \x a — y a \ are equally 
probable and hence (\x a — y a \) = 1 — u a /2. In turn, A = 
1 — (u m ) /2 — \/N/ (2tt) which can be easily confirmed by 
numerical computation. 



2.2 Affecting the inequality 

Blind maximization of the total utility may not be the 
best policy because it can result in large inequalities be- 
tween society members. In an effort to prevent that, the 
matchmaker may adopt the utility function 



u' a (x a ,y a ) = (x k a +y k a ) 1,k 



(k > 0). 



(5) 



where the variants that have one or both utilities negative 
are automatically excluded (they do not have the chance 
to become selected anyway). Choosing k 7S> 1 in Eq. ([5]) 
favors those pairs (x ai y a ) where at least one of the utili- 
ties is high, while k < 1 favors more equal splitting of the 
total utility. The expected value of u' m :— max^ =1 u' a can 
be computed in the same way as in the previous section, 
yielding 



r(i + i)4 1 - 1 /fc 



1/2 



(6) 



Interestingly, (u' m ) scales with N in the same way as (u m ). 

The expected buyer-vendor inequality is now a func- 
tion of k, A(k). Numerical results shown in Fig. [1] confirm 
our initial insight that A(k) grows with k. In particular, 
the value 0.5 achieved for k > 100 corresponds to pick- 
ing the variant which maximizes one of the utilities and 
paying no attention to the other utility. In that case, the 
selected variant has the larger of the two utilities close to 
1 and the other one is 0.5 on average: together we have 
linifc^oo A(k) — 0.5. Further, it can be seen in Fig.[T]that 
k < 1 does not significantly decrease the inequality A(k) 
and thus to create a more social society, one has to use 
a different utility function. For example, for N — 1 000, 
optimizing the outcome of the weaker (which corresponds 
to u a (x a ,y a ) = mhi[x a ,y a }) decreases the inequality by 
29% while reducing the total utility by less than 0.3%. 



2.3 Serving several buyers at once 

When there are several buyers in the market, the match- 
maker can either find the best variant for each buyer sepa- 
rately or she can compromise buyers' needs by finding one 
variant for all. While the former case is identical with our 
analysis above, the latter case is different and requires the 
straightforward generalization of the total utility to the 
form My a + Xi^ a where M is the number of buyers 

and Xi^ a is the utility of variant a for buyer i. For con- 
sistency with our previous formalism, we introduce the 
per-buyer utility 



<(xi 



,XM,a,Va) = Va + 
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Fig. 1. Inequality A(k) as a function of k for N = 1 000. The 
inset focuses on small values of k and the results are averaged 
over 10 000 independent realizations. 
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Fig. 2. Average utility with multiple buyers (results are aver- 
aged over 1 000 independent realizations) . 



where a a is the average buyers' utility of object a, the 
maximal utility we denote as u'^. When utilities 
independent and the number of buyers is large, the central 
limit theorem states that the difference a a — (xi :Ct ) is ap- 
proximately a normally distributed quantity with variance 
proportional to 1/M. It follows that due to the fast decay 
of the normal distribution, the matchmaker cannot find an 
object with the average utility a a differing substantially 
from (xi l0t }- This is confirmed by Fig. [2] where we show 
(uj,) for various values of M. As M increases, fluctua- 
tions of a a gets smaller and (u'^j 1, corresponding to 
a negligible contribution of a a to u'^ (it can be shown that 
to achieve non-negligible max Q a a , the necessary number 
of variants is proportional to e M ). 

Now we see that it's impossible to satisfy several buy- 
ers with one variant. Since production of an individual 
variant for each buyer is often too expensive, it is then 
a natural question how to compromise between the buy- 
ers' satisfaction and the costs of personalized production. 
Within the given framework, one can introduce an addi- 
tional cost which increases with the number of variants 
produced by the vendor — such a cost forces the vendor to 
narrow down the selection. This aspect of buyer-vendor 
interactions is extensively studied in [l|J [2(| where they 
show that based on the compromise described above and 



a few simple additional assumptions, one can reproduce a 
rich variety of market phenomena. 



2.4 Correlated utilities 

So far we assumed that the vendor's and buyer's utilities 
are mutually uncorrelated. While convenient for analytical 
computation, this is not a realistic assumption because in 
general: what is good for the vendor is not good for the 
buyer and vice versa. In other words, one expects x a and 
y a to be negatively correlated. 

To study the influence of correlations we first need to 
find a way how to introduce them into the system. In gen- 
eral it is easy to create correlated quantities by introduc- 
ing a control parameter t € [0, 1] and assuming that both 
quantities have a common part which is proportional to 
t and independent parts which are proportional to 1 — t. 
However, when each part itself is uniformly distributed, 
the resulting distribution depends on the value of t and 
this effect distorts further analysis of the system [l9|] ■ This 
motivates us to switch from uniform distributions to nor- 
mal distributions which preserve their functional form un- 
der addition. We assume that utilities x a and y a are ob- 
tained as 



x a = Vl-tX a + VtC a , y a = Vl-tY a + sViC a (8) 

where X a ,Y a ,C a are drawn from the standard normal 
distribution 7V(0, 1), t E [0,1] is the parameter control- 
ling the correlation strength and the parameter s switches 
between positive (s = 1) and negative (s = — 1) correla- 
tions. Since when adding two normal distributions, indi- 
vidual means and individual variances sum up to give the 
resulting mean and variance respectively, both x a and y a 
have zero mean and unit variance. It is simple to compute 
the Pearson correlation coefficient of x a and y a which is 

To study the system of utilities produced by Eq. {8J, 
we assume the linear total utility given by Eq. |T]). Hence 
u a is normally distributed with zero mean and its vari- 
ance can be shown to be equal to 2(1 + st) := v. We are 
again interested in u m := max^ =1 u a and (u m ). It can be 
shown (see Appendix [5] for details) that (u m ) approxi- 
mately solves the equation 



exp 



4(1 + st) 



N 



1 + st 



(9) 



A comparison of this result with a numerical computa- 
tion of (u m ) (where we randomly generate the utilities 
x a and y a , find the maximal total utility u m and average 
over many realizations) is shown in Fig. [3] As we can see, 
positive correlations amplify the variance of u a and hence 
allow the matchmaker to reach a higher optimal utility. 
Notice that by setting t = in Eq. ([9]), one automatically 
obtains the result for normally distributed uncorrelated 
utilities. 

More insight can be gained if we attempt to find an 
approximate solution of Eq. @. When N is large, the 
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Fig. 3. The dependence of the average maximal utility (u m ) 
on st for N = 1 000 (numerical results and their standard de- 
viations are obtained from 1000 realizations of the model). 
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Fig. 4. The total utility of the selected variant with and 
without the matchmaker. Individual utilities are drawn from 
7V(0, 1), results are averaged over 1 000 realizations. 



factor (u m ) on the left side of Eq. ([9]) is much smaller 
than the exponential term and hence it can be neglected. 
The simplified equation can be solved and gives us the 
approximate result 

(u m ) 2 « 4(1 + st) In [Ny/(1 + st)/n]. (10) 

In contrast to Eq. this time (u m ) grows with N with- 
out bounds. On the other hand, this growth is extremely 
slow: (u m ) is proportional to the square root of In A. For 
example, in the uncorrelated case increasing TV from 1 000 
to 1 000 000 increases (u m ) only by 50%. 

3 Trading without the matchmaker 

Despite all the results obtained so far, one question re- 
mains open: what is the matchmaker's contribution to the 
studied vendor-buyer matchings? This question can be an- 
swered by investigating matchmaker-free methods of the 
variant selection. 



3.1 Vendor proposes 

In [TU they assumed that when the vendor offers a variant, 
the buyer accepts it only if his cost is smaller than the ven- 
dor's cost (instead of maximization of utilities, they stud- 
ied minimization of costs). In our framework this means 
that the buyer accepts variant a only if x a > y a (i.e., 
he wants to profit more than the vendor). The vendor's 
advantage is that he decides which variants to propose — 
obviously, it is optimal to begin with the variant that max- 
imizes y a . This concept, where proposing and accepting 
sides are well defined and distinguished, is similar to the 
classical Gale-Sha pley algorithm known from the stable 
marriage problem [15| . 

Based on the matching described above and using Eq. |T]) 
we can compute the average total utility of the selected 
variant and compare it to the case with the matchmaker. 
Assuming normally distributed utilities and zero correla- 
tions, we studied the system numerically. As can be seen 



in Fig. [?1 the total utility with the matchmaker is consid- 
erably higher. On the other hand, the inequality between 
the fellows decreases from approximately 1.1 (with the 
matchmaker) to approximately 0.4 (when the vendor pro- 
poses). When correlations are present, the difference be- 
tween the two matching methods decreases and becomes 
zero for C xv = ±1. 



3.2 Buyer's search 

As we have seen, the matchmaker optimizes the total util- 
ity at the cost of compromising the utilities of individuals. 
The buyer can avoid being "compromised" by searching 
for the best variant by himself. The drawback is that the 
search is costly (it consumes buyer's time and attention) 
and the corresponding cost has to be subtracted from the 
utility of the eventually selected variant. Since the time 
spent by searching grows linearly with the number of ex- 
amined variants, it is natural to assume the linear cost 
term (3N, where A is the number of examined variants 
and (3 > is the cost per examined variant. The expected 
buyer's utility is 

u s (P, N) = < max x a ) - (3N := (x m ) - (3N. (11) 

1<Q< N 

This form of the buyer's utility was suggested in [ll| where, 
however, the emphasis was on the discussion of the cost 
of information. In today's computerized and networked 
world, searching is easy. Hence to obtain approximate an- 
alytical results, we assume that f3 is much smaller than 
the typical utility value. 

In Eq. (jlip . the term (x m ) grows with A but the cost 
J3N eventually takes over and the total utility us(f3, N) de- 
creases (see Fig. [5ji for an illustration). This behavior is in 
agreement with the classical observation "Good things sa- 
tiate, bad things escalate." by psychologists Coombs and 
Avrunin [2l| . It is now natural to ask, what number of 
examined variants A opt maximizes us(/3,N). We shall do 
that in three distinct cases. In addition to x a drawn from 
the uniform distribution U{— 1,1) and from the normal 
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distribution 7V(0, 1), which were studied before, we con- 
sider also the case when x a is drawn from the power-law 
distribution f(x) = (7 — l)x -7 , x S [1, oo), 7 > 2 (for 
accounts on the importance of power-law distributions in 
complexity management and organization see e.g. 

In Appendix [Bl we find expressions for (x m ) in all three 
cases. Therein, the optimal number of variants -/V op t is 



found in the forms 



Nr. 



uniform: iv opt 
normal: N opt 
power-law: N opt 



(/3/2) 



-1/2 



(/V-M27T/32)) \ 

((3( 7 -i)/r(S)y 1/s , 



(12) 
(13) 
(14) 



o uniform distribution 

^ □ normal distribution 

O ~ o power-law distribution x>' 

4- ^ 




io~ 



where (5 =(7 — 2)/ (7 — 1). Noticeably, in all three cases 
we observe a power-law dependency of ./Vopt on (3. In the 
case of uniformly distributed utilities, (x m ) has an upper 
bound and with N it grows very slowly (that means: there 
is little to be gained by an extensive search). In conse- 
quence, N opt is proportional to and hence it is little 
sensitive to changes of (3. When utilities are normally dis- 
tributed, (x m ) grows with N without bounds — this allows 
-/V^t to reach higher values than in the former case. In 
the power-law case, values of N opt are largest and their 
dependency on (3 is strongest. Moreover, as the power-law 
exponent 7 gets closer to 2, the growth of (x m ) with N 
becomes linear and hence when the growth rate is larger 
than (3, it is optimal for the buyer to search "forever" (the 
exponent 8 diverges). 

To review the accuracy of the presented results, in 
Fig. [5] we compare them with numerical simulations. As 
can be seen, good agreement is achieved in all three cases. 
Let us conclude with the remark that in the power-law 
case, the average, the mode, and the median of x m may 
differ significantly. In consequence, it is an important ques- 
tion whether the buyer should rely on (x m ) which is strongly 
affected by rare extreme events. 



4 Conclusion 

In this paper, we studied several simplified models of inter- 
actions between buyers and a vendor. Assuming that both 
the vendor and the buyers can attribute certain utility val- 
ues to each available variant, the decision which variant 
is to be delivered can be formulated as a mathematical 
optimization problem. 

Firstly we assumed that there is a matchmaker who 
can fairly select the optimal variant according to what is 
best for all participants of the contract. A plausible cri- 
terion for the matchmaker is to select the variant with 
the maximal total utility. Our results show that when the 
utilities are uniformly distributed, the matchmaker can 
achieve the total utility close to its upper bound even 
when the number of available variants is small. In other 
words, little choice is enough in this case. On the other 
hand, when the total utility is maximized, the difference 
between utilities of the vendor and the buyer may be too 
large to consider it to be the optimal choice. When the 
mere summation of individual utilities is replaced by a 
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Fig. 5. The buyer's search for three distinct probabilistic dis- 
tributions of utilities — comparison of numerical results (shown 
with symbols) and analytical results (shown with lines) for 
«s(0.01, N) (a) and N opt (b); 7 = 4 in both figures. 



more refined expression, the inequality between the in- 
volved parties may be decreased. In particular, we found 
that when the matchmaker maximizes the smaller of the 
two individual utilities, the inequality decreases signifi- 
cantly while the total utility is almost unchanged. 

When the matchmaker tries to find one variant for 
several buyers, the situation turns out to be far less fa- 
vorable: the number of variants needed to approach the 
upper bound of the total utility grows exponentially with 
the number of buyers. As a result, production of multiple 
variants is advisable — this topic was extensively studied 
in our previous works [T^ . [20 |. Assuming normally dis- 
tributed utilities, we studied the influence of correlations 
on the system behavior. While analytically more demand- 
ing, this generalization is important because it makes the 
system more realistic. Our results confirm that the opti- 
mal total utility depends strongly on the correlation of 
utilities. 

Secondly we studied the variant selection without a 
matchmaker. In this case, the space of possible means of 
variant selection is vast and hence we focused on two par- 
ticular situations. In the first one, the vendor offers and 
the buyer passively decides whether to accept the offered 
variant or not. While in comparison with the matchmaker- 
mediated outcome, the total utility is slightly lower (see 
Fig. the inequality is decreased substantially. In the 
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second selection method, the buyer chooses the variant 
by himself but he also has to pay the cost for examining 
the variants. Here we considered three different distribu- 
tions of utilities and shown that when the distribution has 
a power-law tail, for the buyer it may be optimal to in- 
spect a huge number of variants. In the extreme case of a 
power-law distribution with the exponent lower or equal 
than 2, the buyer does best by searching forever. Admit- 
tedly, these results are influenced by the linear growth of 
searching cost with the number of examined variants N 
assumed by Eq. () 11 [1 . In contrast, in psychology it is well 
known that refusing the second best variant or having too 
many options may be frustrating for people [231 ] . These ef- 
fects could be included by adding an additional cost term 
depending on the utility of the second best variant, or 
simply by a part of the searching cost proportional to a 
higher power of N . 



A Extreme statistics for normally distributed 
variables 



In this appendix we show how (u m ) can be approximated 
when the number of variants, N, is large. Assuming u a £ 
u m := max*_j u a has the distribution 



f(u m ) = 



Nc -u 2 m /2v 



2 /2v 



N-l 



-du 



(15) 



The error function erf [x] := -^U J* exp[— t 2 ] dt allows us to 

write the integral in Eq. (|15p as \ — erf [u m / \/2v]/2. When 
N is large, u 2 n /(2v) ^> 1 and hence we can use the asymp- 
totic expansion erf [x] ps 1 — exp[— x 2 ]/(x^/tt) (taken from 
matliworld.wolfram.com, the next contributing term is 
proportional to exp[— x 2 ]/x 3 ). Again, we use the approx- 
imation 1 — y ps exp[— y] (which is valid for y -C 1) to 
obtain 



: exp 



2v 



Unfortunately, the integral corresponding to (u m ) cannot 
be solved. Since f(u m ) does not have heavy tails, a rea- 
sonably precise result can be obtained by approximating 
(u m ) with u m maximizing f{u m ), yielding 



(u m ) = (N - 1)(1 + v/(u m ) 2 )J — exp[-(M m ) 2 /24 

When N is large, vj (u m ) 2 <C 1 and hence we neglect this 
term on the right side. In addition, N — 1 « N and we get 



(u m ) exp [(u m f/2v] = N ^/2tt (16) 
which after substituting v = 2(1 + st) gives Eq. 



B Analysis of the buyer's search 

We begin with x a G U(— 1, 1). Then x m has the distribu- 
tion g(x m ) — y [l — (1 — x m )/2\ and in consequence 
(x m ) = 1 — 2/{N+ 1). Maximizing us{(3, N) with respect 
to N we get A opt = t/2/P - 1. 

Now let's consider x a e jV(0,1). After substituting 
v = l in Eq. P^|) . the expected utility ug(P,N) can be 
maximized using implicit derivative, yielding the condi- 
tion (3N{x m ) = 1. In this equation, "fast" and "slow" 
terms {N and (x m ) respectively) are mixed together and 
hence an approximate solution can be found by the fol- 
lowing iterative procedure. Approximating (x m ) = 1 gives 
the rough estimate Nq = 1/(3. Together with Eq. (TIB"]) , this 
value leads to the improved estimate (x m ) ps ■J— ln(27T/9 2 ) 
and in turn we get 



A opt « (/V-M2tt/3 2 )) \ 

an improved estimate of N opt . 

Finally let's study the case where utilities x a are con- 
strained to the range [1, oo) and follow the power-law dis- 
tribution f(x a ) — (7— l)i' _7 , 7 > 2. Then the cumulative 
distribution is P(x a ) — x 1 ^ 1 and the distribution of x m 
is g(x m ) = AT(7— l)x r (l — a; 1-7 )^ -1 . After approximat- 
ing 1 — x 1 ^ 1 ps exp[— x^7 7 ] (the relevant values of x m are 
large) and replacing the lower integration bound with 
we obtain (x m ) ps AtVCt-i)/-^] where s = ( 7 _ 2 )/( 7 _i) ; 
this result is well defined for all 7 > 2. When only a 
fraction r of all utilities follows the power law, the re- 
sult generalizes to (x m ) ps (Ar) 1 ^ 7-1 ) r[5]. Maximization 
of the expected utility us(f3,N) is straightforward and 

yields A opt ps [^(7 - l)/r[<5]] . For a different analysis 
of the largest value of a power-law distributed sample and 
a broader discussion of power laws see [24| . 
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